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Abstract 

A non-perturbative calculation of the gyrotropic pressures associated with large-scale 
mirror modes is performed, taking into account a finite, possibly anisotropic electron 
temperature. In the small-amplitude limit, this leads to an extension of an asymptotic 
model previously derived for cold electrons. A model equation for the profile of subcritical 
finite-amplitude large-scale structures is also presented. 
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PACS: 52.35.Py, 52.25.Xz, 94.30.cj, 94.05.-a 



I. INTRODUCTION 

Pressure-balanced magnetic structures in the form of strong magnetic enhance- 
ments (humps) and depressions (holes) that are quasi-stationary in the plasma 
frame, with no or little change in the magnetic field direction, are commonly ob- 
served in regions of the solar wind and of planetary magnetosheaths with relatively 
large /3 and a dominant (generally ion) temperature in the transverse direction ( 
see, for instance, and references therein). The origin of these structures is 

still not fully understood, but they are usually viewed as nonlinearly saturated 
states of the mirror instability (MI) discovered by Vedenov and Sagdeev j^. It is 
a kinetic instability whose growth rate was first obtained under the assumption 
of cold electrons, a regime where the contributions of the parallel electric field E\\ 
can be neglected. However, in realistic space plasmas, the electron temperature 
can hardly be ignored The linear theory retaining the electron temperature 
and its possible anisotropy, in the quasi-hydro dynamic limit (which neglects finite 
Larmor radius corrections), was developed in the case of bi-Maxwellian distribu- 
tion functions by several authors [sl-jol. A general estimate of the growth rate 
under the sole condition that it is small compared with the ion gyrofrequency (a 
condition reflecting close vicinity to threshold) is presented in [l^ . The instability 
then develops in quasi-perpendicular directions, making the parallel magnetic per- 
turbation dominant. This analysis includes in particular regimes with a significant 
electron temperature anisotropy for which the instability extends beyond the ion 
Larmor radius. In the limit where the instability is limited to scales large compared 
with the ion Larmor radius, only the leading order contribution in terms of the 
small parameter 'j/ {\k\zV\\i) is to be retained in estimating Landau damping, and 
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the growth rate is given by 
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where 



Tu (^11 + 2^1,(^2 ^1) . 1 



2^||(l + ^x)(^|| + l) A 
measures the distance to threshold and 



-^-^ (2) 



E 



2 + 9^{A + 9^) + 9l 



F = ^11- I - 1 + - 
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Here, T±a and T||q are the perpendicular and parallel (relative to the ambient 
magnetic field Bq taken in the z direction) temperatures of the species a {a = 
i for ions and a = e for electrons ), 9± = T^e/T±i, 9\\ = T\\e/T\\i and = 
+ /3±e with f3^a = 8np±a/Bl where p±a is the perpendicular thermal pressure 
(similar definition for Furthermore, the parallel thermal velocity is defined 
as v\\a = \/'2T\\a/ma, and = (2Tij/mp)^/^/r2j denotes the ion Larmor radius 



{Qi = eBo/niiC is the ion gyrofrequency) . 

The growth rate given by Eq. ([1]) has the same structure as in the cold electron 
regime considered in ll| in the case of bi-Maxwellian ions and generalized in 12 1 
and js] to an arbitrary distribution function. The first term within the curly 
jrackets provides the threshold condition which coincides with that given in jsl- 
Uj. The second one refiects the magnetic field line elasticity and the third one 
(where F depends on the electron temperatures due to the coupling between the 
species induced by the parallel electric field which is relevant for hot electrons) 
provides the arrest of the instability at small scales by finite Larmor radius (FLR) 
effects. 
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An aim of this letter is to extend to hot electrons the weakly nonlinear analysis 
previously developed for cold electrons Q- Since in this asymptotics, FLR 
contributions appear only at the linear level, the idea is to use the drift kinetic 
formalism to calculate the nonlinear terms. We show that the equation governing 
the evolution of weakly nonlinear mirror modes has the same form as in the case of 
cold electrons. In particular, the sign of the nonlinear coupling coefficient that pre- 
scribes the shape of mirror structures, is not changed. This equation is of gradient 
type equation with a free energy (or a Lyapunov functional) which is unbounded 
from below. This leads to finite-time blowing-up solutions 15|, associated with 



the existence of a subcritical bifurcation 
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14| . To describe subcritical station- 



ary mirror structures in the strongly nonlinear regime, we present an anisotropic 
MHD model where the perpendicular and parallel pressures are determined from 
the drift kinetic equations in the adiabatic approximation, in the form of prescribed 
functions of the magnetic field amplitude. 



II. BASIC EQUATIONS 

A main condition characterizing mirror modes, at least near threshold, is pro- 
vided by the force balance equation 

^( B\ r 471, ,n(B-V)B 

-V(px + — )+ 1 + ^{P±-P\\' 



67r 

+B(B-V)(^^)-V-n = 0, (3) 

where the pressure tensor, viewed as the the sum of the contributions of the var- 
ious species, has been written as the sum of a gyrotropic part characterized by 
the parallel (py = J2aP\\a) and perpendicular {p± = J2aP±a) pressures, and of a 
gyroviscous contribution 11 originating from the sole ion FLR effects when concen- 
trating on scales large compared with the electron Larmor radius. As mentioned 
above, FLR effects arising only at the linear level with respect to the amplitude of 
the perturbations, the other linear and nonlinear contributions can be evaluated 
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from the drift kinetic equation for each particle species 

^^"+^ b-v/„ + f-/xb-vi? + ^i5;iil|^ = o (4) 

L rrir, J af II 



dt 



We ignore the transverse electric drift which is subdominant for mirror modes. 
In this approximation, both ions and electrons move in the direction of the magnetic 
field (defined by the unit vector b = B/i?) under the effect of the magnetic force 
II h ■ VB and the parallel electric field E\\ = —h ■ V0 where the magnetic moment 
H = v'j_/{2B) is an adiabatic invariant which plays the role of a parameter in Eq. 
(jlj). Here is the electric potential. The quasi- neutrality condition Ue = rii = n, 
where Ua = B J fadiidv\\dip = J fad^v, is used to close the system and eliminate 

^11- 

In this framework where FLR effects are neglected, the gyrotropic pressures are 
given by Pa\\ = rriafvffad^v = m^^B ^ vpadlidv\\dip, and p„_L = f vlfad^v = 

TTlaB'^ J flfadfidV\\dip. 

The asymptotic equation governing the mirror dynamics near threshold is ob- 
tained by expanding Eqs. ([3]), (jl]) and the quasi-neutrality condition, with the 
pressure tensor elements for each species computed near a bi-Maxwellian equilib- 
rium state characterized by the temperatures T±a and T\\a- 

III. LINEAR INSTABILITY 

Before turning to the nonlinear regime, we briefly review the derivation of the 
MI linear growth rate in the simplifled framework provided by the drift kinetic 
approximation which is only valid at scales large enough for FLR effects to be 
subdominant. 

Linearizing Eq. ([3]) about the background fleld Bq and equilibrium pressures 
p^l^ and p^^\ and considering perturbations B and p^f^ oc e~*^*+*^ '", we get 



(1) BqB, _ K ^ /3± - \ BgB, 



Here, p^f' has to be calculated from the linearized drift kinetic equation 



dt 



dz 



+ 



- /^l^ + —^11 -IT- 
oz nia OV\\ 



0, 
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where we assume each /^°) to be a bi-Maxwellian distribution function 



f (0) = A 



exp 



V 



2 



T, 



(7) 



with Aq, = noma/ {271 y/TTV\\aT±a)- 

Equation ([6]) is solved in Fourier representation, as 

° u — kzV\\ dv\\ 



(8) 



The neutrality condition allows one to express the potential in terms of Bz. In- 
deed, assuming ( = ^/^llJJ/{\kz\v\\i) ^ 1 (so that the contribution from the Landau 
pole is small). 



fi^^dvzdfidip 



no 



Bon, 



nn Bz 



(9) 



Similarly, neglecting the electron Landau resonance contribution because of the 
small mass ratio. 



no 



BoT\\e 



rr Bz 
-DO 



Consequently, 



T 



±i 



Oil 



Bz 
Bo 



(10) 
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1 + ^11 

We thus recover that for mirror modes, the parallel electric field vanishes when the 
electrons are cold (6*^ = 6\\ = 0). Interestingly, when 9± = 9\\, only the Landau 
pole contributes to 0. 

It is now necessary to evaluate 

.Bz 



p± = ^^p± + B', J2^aJ ^ifi'^d^Idv^^dcp. 



(12) 
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Using 



u — kzV 

kyVn 



CO — kzV\\ -Do 



-Do 

no 



we get 



+ r + 



Bo' 



(13) 



■ ' ■ Til, 2(1 + 

Substituting this expression into the hnearized force balance equation yields the 
linear instability growth rate given by Eq. ([1]), up to the FLR term which is not 
captured by the drift kinetic approximation. Note that the growth rate given by 
Eq. ([T]) is consistent with the apphcability condition 'y/\kz\ ^ v\\i near threshold 
(r ^ 1), as kz and {kz/k±Y scale like F, while 7 like T^. 



IV. GENERAL PRESSURE ESTIMATES 

As demonstrated in [isl [l^ , the scalings resulting from the linear theory near 
threshold imply an adiabaticity condition to leading order. It is thus enough to 
consider the stationary kinetic equation 

dfa 



t;||b-V/,-(b- V) 



TTLr, 



dv\\ 



0. 



(14) 



It turns out that Eq. f[T4j) is exactly solvable, the general solution being an arbitrary 

function /q, = (7a(/i, Wa) of the particle energy = f + jj,B H —cj), and of /x. 

To find the function ^'^(/i, Wa), we use the adiabaticity argument which means that, 
to leading order, function of /i and Wa retains its form during the evolution. 

Therefore, the function (y'a(/i, Wa) is found by matching with the initial distribution 



function f^^ given by Eq. ([7]) which corresponds to 
We get 



and Wa 



+ l^Bo 



ga{fJ',Wa) = Aaexp 



nBoiria 



Aa exp 



2Wa 



liBoniai^ 



||a 



(15) 
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Thus, QaifJ'jWa) is a Boltzmann distribution function with respect to Wa but, at 
fixed Wa, it displays an exponential growth relatively to fi if T±a > T\\a- This effect 
can however be compensated by the dependence of Wa in /i. This means that only 
a fraction of the phase space (/x, Wa) is accessible, a property possibly related with 
the existence of trapped and untrapped particles. 

Note that expanding Eq. f|T5|) relatively to B^/ Bq and jT^i reproduces the 
first order contribution to the distribution function given by Eq. (|8]) with w = 0, 
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14j | in the case of cold electrons. 



and also the second order correction found in 
It should be emphasized that Eq. ( JT5l) only assumes adiabaticity and remains valid 
for finite perturbations. 

The function can also be rewritten in terms of fy, and as 

,2 , 



ga = Aa exp 



rriaVii 
2T1 



n. 



X 



2 .T 



exp 



Ti 



II" 



2T±a T\\a B 

which can be viewed as the bi-Maxwellian distribution function with the renormal- 
ized transverse temperature 

-1 



T 



T, 



T, 



Til, 



_ Bq / Tj_a ^ 



||a 



(16) 



Note the Boltzmann factor exp —[ea(f)/T\\a] in the expression of g^- For cold elec- 
trons, the ion distribution function was obtained in 16|] by assuming that it remains 
bi-Maxwellian, and owing to the invariance of the kinetic energy and of the mag- 
netic moment. This estimate, obtained by neglecting both time dependency (and 
consequently Landau resonance) and finite Larmor radius corrections, reproduces 
the closure condition given in 17 1. 
After rewriting Eq. f lTSj) in the form 

1 „.2 



ga = Aa exp 



Til 



flBQlTLa 
T±a 



(1 + 



T\ a B — Bn 



Ti, 



Bn 



(17) 
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the quasi-neutrality condition gives 

T±i B — Bq 



1 



1 + 



T\\i Bq ^ 

T±e B — Bp 
T\\e Bo 



( 



exp 



ecj. 



or 



e0=(rri + T|r;)-ix 



log 



1 + 



T \ f. B — Bn 



fin 



1 + 



T \ i B — B(\ 



-1 



2ih 



fin 
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Interestingly, the electron density (and thus also that of the ions) 



Hp 



B A , Tj_e B — Bq 



fio 



Til, 



fin 





60 


1 exp 


Til 



(19) 



has the usual Boltzmann factor exp e0/T||e and also an algebraic pref actor de- 
pending on the magnetic field fi. In the case of isotropic electron tempera- 
ture (T^e = T||e = Te), the electron density has the usual Boltzmann form 

He = no exp [e(f)/Te]. 

Equation (fT8|) shows that the potential vanishes in two cases: for cold electrons 
and also when electron and ion temperature anisotropies and (with = 
T\ nlT ^n) are equal, a case considered in the linear theory of the mirror instability 
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In order to evaluate explicitly the perpendicular pressure for each species 



rriaB^ 



fi, 



fiQadfidviidLp 

T_La B — Bq 

Tiu fio 



exp 



Ti, 



where ecj) is given by Eq. (jTSj) . it is convenient to introduce the functions 



S I Ju) 



1 + aiu 
1 + u 

1 + tteU 
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1 + OgM/ 
^ ^1 + tteUY" 



1 + aiU 



(20) 
(21) 



with the notations u = (B — Bq)/ Bq and Cq, = / {Ti^i" + ^||e^)- The two latter 
functions transform one into the other by exchanging the subscripts i and e. The 
ion and electron perpendicular pressures are then written as p^^ = noT^aSLaiu). 
In the special case of cold electrons, 

B^ , rj_j B — Bo 

' ||i ^0 



P. = n„r,„(i + ^— ^) , (22) 



which is algebraic relatively to B. From this expression as well as from the gen- 
eral formula for p± = p±i + p±e given by Eqs. (l20l) and ( |2T]) it follows that the 
perpendicular and magnetic pressures are anticorrelated. When B increases (de- 
creases), the ratio of the perpendicular to the magnetic pressure, i.e. the local 
/3±, decreases (increases), which corresponds to a reduction (an increase) of the 
distance to threshold. This implies that the instability cannot saturate at small 
amplitudes. 

Similarly, for the parallel pressure, we have 

^ -B / T_i_aB-Bo\-i / ea(j)\ . . 

that rewrites p\\a = '"'0^||a'S'||Q,(M) with 

S,(u) = (^) (i±^y (24) 

VI + aiuj Vl + aeUj 

S„iu)=(^:±^)C-±^y\ (25) 

V 1 + tteUj V 1 + ttiUj 

V. THE WEAKLY NONLINEAR REGIME 

As it follows from Eq. in the linear regime near threshold, the fluctuations 
of perpendicular and magnetic pressures almost compensate each other. In the 
weakly nonlinear regime, the second order correction to the total (perpendicular 
plus magnetic) pressure is thus relevant and leads to a local shift of F. To find this 
correction, we consider the expansions of the perpendicular pressures of the ions 
and electrons in the u variable. Because of the symmetry between the functions 
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S±i{u) and S±e{u), it is enough to consider the expansion 



S±i{u) = 1 + u{2 ~ 2ai - Ci{ae - Oj)) 



-4aj + 3a^ + -(?i{ae 



+2aiCi(ae — ttj) + 1 + O (u 



2ci(ae - tti) 



As a result, the second order contributions to the perpendicular ion pressure is 
given by 



3a,- - 4aj + 1 + Qfoe - aA 



X [2('^i + l)(«e - flj) - 2 + 3a, 



with an analogous formula for the perpendicular electron pressure, obtained by 
exchanging the i and e indices. Furthermore, the threshold condition rewrites 



+ no\T±i [2 - 2ai - q (og - a^)] 
47r ^ 

+Tj_e [2 - 2ae + Ce (tte - fli)] } = 0. 



(26) 



The quadratic contributions to the pressure balance (j3]), originating from p\/ + 
P?± + {.B - Bof /(Stt), are collected in a term A {^^Y with 

A = no{T±i{3a^^ - 4ai + 1 

hi + Ci){ae - ttj) - 2 + 3ai 



-\-Ci[CLp do 



+T±e[3al - 4ae + 1 + Ce{ae - tti 



X 



^2 

;i + Ce)(ae - a^) + 2 - 3aJ)| + 



(27) 



The value A^ of A at threshold is obtained by expressing Bq/8tt by means of Eq. 
( I26l) . which gives 



= nQ^T±i 3a^ — 4aj + 1 
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+Ci{ae - + Ci)(ae - a^) - 2 + 3a, 



--(2 - 2ai - Ci{ae - 



Sal — 4ae + 1 + Ce(ae — flj) 



X 



(-(1 + Ce)(ae - ttj) + 2 - 3ae] 
-i(2 - 2ae + Ce(ae - ttj)) }. 

After some algebra, defining Ac = ^c/{noT^i), one gets 

2 



A 



c 



Til,; 



^3 

3 _|_ 3^ _ _ 



2^11 (l + ^| 



, 2 



(^± + ^||) + 2^11(1 + 



(28) 



Proceeding as in 13|, retaining the contribution of the above quadratic terms 
to the pressure balance, leads one to supplement a nonlinear contribution to Eq. 
([U) that becomes 



du 
'dt 



2 T\\i vi\i 
V^T^^ D 



/c,{r 



A. 



, 3 /Tj_i X 1 + F 2 , 

+4(7i:-^)rMi'-^-"-2(i+»,) 



u 



(29) 



Here the integral operator /C^ reduces in Fourier representation to \kz\ and x 
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. Furthermore, within the present approximation, u coincides with B^/ Bq. 
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Equation fl29|) extends the result of 
latter case, this equation is a gradient type equation 

5F 



2l| valid for cold electrons. As in the 



du 
dt 



Su ' 



for which the free energy (written in dimensionless variables) 
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is unbounded from below due to the integral / XcU^dr. This leads to a blow-up 
behavior, associated with a subcritical bifurcation [isl, [l^. Saturation at large 
values of the amplitude and formation of stationary structures requires additional 
nonlinear effects such as the influence of resonant particles on the nonlinear cou- 



pling [22|. Equation (129|) that does not include saturation processes is not suitable 



to address the question of the reduction of the temperature anisotropy by the devel- 
opment of the mirror instability mentioned in jsj. This effect is reproduced by the 



quasi-linear theory 19|, and was also studied in the context of the so-called FLR- 
Landau fluid model that, like the present asymptotics, retains a linear description 
of the Landau resonance and of FLR effects, but includes all the hydrodynamic 
nonlinearities and does not a priori prescribe a pressure balance condition. It was 
observed in this case that during the saturation phase, the mean temperatures 



2Q|. 



rapidly evolve in a way as to reduce the distance to threshold 

As demonstrated in [l^, the sign of the nonlinear coupling Ac defines the 
type of the mirror structures, namely holes (Ac > 0) or humps (Ac < 0), near 
;hreshold. This sign is strongly dependent on the equilibium distribution function 



23( 1 . It is nevertheless of interest to consider the case where both ions and electrons 
have a bi-Maxwellian distribution function. It turns out that the sign of Ac can 
then be determined analytically in a few special cases. 



(i) Limit 6\\ < 9±: 



Ac e^, fT±e 3\ , , 

^ ' - - > 0. (30) 



noTj_iai 9\\ \T\ 
(a) Equal anisotropics (9± = 9\\) 

Ac = no{T^i + T^e) (3a2-4a + l 

-no{T^i + Txc) (1 - a) = 3a-^ > 0. 

OTT 

(Hi) Isotropic electron temperature: The coefficient Ac can be rewritten in the form 
Ac = no{ai - l){T±i({3ai - 1) 
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+Ci [-(1 + Ci) (aj - 1) + 2 - 3ajJ ^ 
+TeCe[i(l + Ce) (ai-l) + l]} + 



) 




Furthermore, at threshold, 




(31) 



Hence, we simultaneously have two inequalities > 1 and TieCe > T±i{ci — 2). 
Therefore, 



which is positive, because 1 — q = Ce = (1 + 6*11)"^ > and Oj > 1. 
(iv) More general conditions: A numerical approach was used in this case. Figure 
1 displays, for typical values of the parameters (taken here as 6'_l = 1, = 1.1 
and = 10), the distance to threshold F (dashed line) given by Eq. ([2]) and the 
non-dimensional nonlinear coupling coefficient A = A/{noT±i) (solid line), where 
A is given by Eq. (127|) . as a function of 6\\. This graph is typical of the general 
behavior of these functions and shows that they are both decreasing as 6\\ increases, 
with A possibly reaching negative values, but only below threshold. In order to 
show that the value Ac, given by Eq. (128|) . of A at threshold is positive in a wider 
range of parameters, we display in Fig. 2, as a function of P±i for 6± = 0.2 
(solid line), 6± = 1 (dotted line) and 9± = 5 (dashed line), the quantity min (Ac) 




+T,Ce[-{l + Ce){ai-l) + l\} 

+ ^no{ai - 1) [T^i (2 - d) + T^eCe] 
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FIG. 1: Fig. 1. Variation with 6\\ of the distance to threshold F given by Eq. ([2]) (dashed 
line) and of the normalized nonlinear coupling coefficient A (solid line) evaluated from 
Eq. dSID for 6l_L = 1 , ai = 1.1 and I3±_i = 10. 

1 0.00 \ i 




FIG. 2: Fig. 2. Variation with f3±_i of the minimum min (Ac) of the normalized nonlinear 
coupling coefficient taken in an interval of values of ap between and api(/3j_i), defined 
such that the threshold is obtained for a value of 6\\ equal to 100, for 9±_ = 0.2 (solid 
line), = 1 (dotted line) and 9± = b (dashed line). 

obtained after minimizing Ac in an interval of values of ap between and api(/3±i). 
The latter quantity is arbitrarily defined such that the threshold is obtained for a 
value of equal to 100. This graph shows that min(Ac) varies little with ^_|_ but 
is very sensitive to P±i. As the latter parameter is increased, min (Ac) decreases 
but remains always positive. Although this numerical observation is not a rigorous 
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proof, it convincingly shows that A > in the parameter range of physical interest. 



VI. STATIONARY NONLINEAR STRUCTURES 



Substituting the explicit expressions of the gyrotropic pressures in terms of 
the magnetic field amplitude given in the Section IV, within the equation for the 
balance of forces 



B 



An 



(B ■ V)B 

47r 



-B(B- V)( 



0, 



(32) 



leads to a closed system that seems overdetermined due to the divergenceless con- 
dition V ■ B = 0. In fact, it can be checked, after some algebra using the explicit 
expressions ( 120|2T]) and ( 12^11251) . that the projection of Eq. ( l32l) on the magnetic 
field vanishes identically, thus reducing the system to three equations for three 
unknowns. These equations can be useful for finding, possibly numerically, sta- 
tionary profiles of three-dimensional finite-amplitude stationary mirror structures. 
Note that Eq. (!32|) differs from the Grad-Shafranov equation 2J, |25|] in that the 
parallel and perpendicular pressures are here prescribed functions of the magnetic 
field amplitude. A main issue concerns the existence of stable subcritical solutions, 
a question that is beyond the scope of this letter and will be addressed in forthcom- 



ing works. Such structures are reported by satellite observations 
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also expected from the subcritical character of the mirror instability 



27| and are 



14. Equilib- 



17|, where they lead to 



rium solutions were computed in one-space dimension in 
discontinuous profiles. Their regularization would require that FLR corrections be 
retained. These additional contributions are known from the linear kinetic theory 
but their extension to the finite-amplitude case remains a challenging problem. 
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